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We introduce a method that can orthogonalize any pure continuous variable quantum state, i.e. 
generate a state \ip±) from \ip} where {^p\tp±) = 0, which does not require significant a priori knowl- 
edge of the input state. We illustrate how to achieve orthogonalization using the Jaynes-Cummings 
or beam-splitter interaction, which permits realization in a number of systems. Furthermore, we 
demonstrate how to orthogonalize the motional state of a mechanical oscillator in a cavity optome- 
chanics context by developing a set of coherent phonon level operations. As the mechanical oscillator 
is a stationary system such operations can be performed at multiple times, providing considerable 
versatility for quantum state engineering applications. Utilizing this, we additionally introduce a 
method how to transform any known pure state into any desired target state. 
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A qubit basis formed by a pair of orthogonal quan- 
tum states is central to quantum information processing. 
Currently there is considerable effort towards implement- 
ing quantum information processing with two-level sys- 
tems. For such systems, an intriguing and fundamental 
fact is that quantum mechanics prohibits the construc- 
tion of a universal-NOT gate that would produce an or- 
thogonal qubit from any input qubit [1]. This quantum 
mechanical property is closely related to the quantum 
no-cloning theorem [5], however, faithful cloning can be 
achieved probabilistically provided that the set of input 
states is linearly independent J3 . Similarly, using such 
an input set of states, it is possible to construct a proba- 
bilistic NOT operation for qubits [1]. A qubit basis may, 
however, also be formed using two orthogonal continuous 
variable states. Thus far, efforts to construct such a basis 
have mainly concentrated on using a superposition of co- 
herent states [S]. Also, recently a qubit basis was realized 
using photon subtraction from squeezed vacuum ^ . 

In this Letter, we introduce a method for quantum 
state orthogonalization for continuous variable quantum 
systems. Notably, the method only requires knowing 
the angle i9 made by the state's mean amplitude (6) — 
1(6)1 e"', where b is the annihilation operator, and hence 
the scheme is magnitude independent. Furthermore, our 
method is readily extended to generate an arbitrary su- 
perposition of the initial state and an orthogonal coun- 
terpart to allow the encoding of quantum information. 
The orthogonalizer T|_ ex be~'''^ + b^e^'^ is formed by a lin- 
ear superposition of the bosonic annihilation and creation 
operators and generates a state orthogonal to any pure 
state \ip) i.e. (V'| T_i_ |V') = when (j) = d + 7r/2. Thus, Tj_ 
is a quadrature operator that is perpendicular to t? [7]- 

The orthogonalizer can be realized with interactions 
that are available in many physical systems, e.g., to re- 
alize Tj_ in cavity-quantum-electrodynamics [5J [S], one 
prepares an input qubit in the state A\g) + B |e) which 
then weakly interacts via the Jaynes-Cummings Hamil- 
tonian H/h = —in{b(T^ — b^a^), where fl is the coupling 



rate and ct-i-.- are the raising and lowering operators. A 
controUably weighted superposition of addition and sub- 
traction is achieved by projective measurement of the 
qubit onto B* \g) —A* |e). The measurement operator is 
thcnTQED = i{g\B-{e\A)il-nT{ba+-b^a.mA\g) + 
B |e)) = nriA'^b+B'^b'^), seeFig.l (a). With this interac- 
tion, optical [S] or microwave [3] fields in a cavity, or the 
motional state of trapped ions [10], can be orthogonal- 
ized by appropriately setting A and B. Similarly, a pure 
state of a traveling optical field can be orthogonalized by 
interaction on a beam-splitter and then measurement of 
an optical qubit comprising a superposition of zero and 
one photons JJJ, see Fig.l (b). As these interactions are 
common throughout quantum optics, adaptations of this 
orthogonalization protocol to other physical systems can 
be readily achieved. Moreover, a different scheme to per- 
form a superposition of photon subtraction and addition 
was recently proposed fT^, which could also be used to 
realize state orthogonalization. 

The tools we introduce for orthogonalization can also 
be utilized for quantum state engineering applications. 
Currently, single-quanta-manipulation techniques per- 
formed on traveling light fields 13J have prepared su- 
perposition states via photon subtraction 1 14', observed 
the bosonic commutation relation 15J, and engineered ar- 
bitrary quantum states up to the two-photon level jl6) . 
Much progress has also been made for arbitrary quantum 
state preparation of the motion of trapped ions and mi- 
crowave field states [T7j . As mechanical elements are now 
also considered for quantum applications, experimen- 
tal tools are required for the coherent manipulation of 
phononic modes. Examples of progress in this direction 
are the observation of the ground state of motion [T5U20J. 
steps towards single-phonon manipulation by coupling to 
a superconducting phase qubit [IH] , strong coupling ^T\ , 
and mechanical mode thermometry via sideband asym- 
metry ^224. Also, recently the lattice vibrations of two 
diamonds were entangled by coherently distributing one 
quanta across the two vibrational modes ^5] . 



Coherent phonon manipulation. — In this section we 
demonstrate how to perform an arbitrary coherent super- 
position of phonon subtraction, addition and the identity 
operation to a mechanical oscillator using cavity optome- 
chanics. The prototypical optomechanical system is a 
Fabry-Perot cavity where one of the mirrors is sufficiently 
compliant that the reflection of light can modify the mir- 
ror momentum via radiation-pressure. Concurrently, the 
motion of the moving mirror modulates the optical phase 
and generates sidebands. To realize phonon subtraction 
(addition) one can optically drive an optomechanical cav- 
ity at the red (blue) sideband and then perform single 
photon detection on the field scattered onto cavity res- 
onance. Provided that the sidebands are well resolved 
and the optical phase shifts are small allowing lineariza- 
tion, the red-detuned drive gives rise to a beam-splitter 
interaction and the blue-detuned drive gives rise to a two- 
mode-squeezing interaction. This linearization procedure 
was discussed, for optomcchanics, in Ref. [24j, where 
quantum state transfer between light and mechanics was 
proposed. Drive on the blue sideband has also been con- 
sidered for continuous-variable teleportation from light 
to the mechanics |25) . Some other applications utilizing 
these sidebands are reviewed in Refs. ^ . 

Our proposed setup for coherent phonon control uses 
two orthogonally polarized optical fields to interact with 
the mechanical resonator, see Fig. 1(c). We consider a 
pulsed protocol where the conditional mechanical state 
following the pulsed interaction and measurement is de- 
termined. The optomechanical Hamiltonian [27 for the 
two independent optical modes in the optical rotating 
frame at the drive frequencies is 
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where Hd/h = J2i V2i^Ni{£*ai +£ial) is the drive term, 
the subscripts label the two orthogonally polarized modes 
i€{h,v}, and a (b) is the cavity (mechanical) annihi- 
lation operator, (wm, mechanical angular frequency; 
A, optical detuning; go, optomechanical coupling rate; 
K, cavity amplitude decay rate; N, photon number per 
pulse; £, drive amplitude, where /dt|£| =1.) Neglect- 
ing mechanical damping and input noise, as the inter- 
action time can be made shorter than the decoherence 
time scale, we compute the dynamics in a similar man- 
ner to Ref. [21] • The mechanical evolution is computed 
via the Hamiltonian and the cavity field is computed via 
the Langevin equation Oi = —ia.i [A^ — gQ{b + b^)] + 



v2K(ai, 



fNiEi) — KOi, where Oin is the optical input 



noise. We enter a displaced frame to follow the mean of 
the operators, i.e. Oi -^ ^/NiUi + Oi and b ^ (H + b. Pro- 
vided that the intracavity intensity varies much slower 
than the mechanical frequency the mechanical mean am- 
plitude is /3 ~ ^'l2i^i\'^i\ ■ This displacement due 
to the optical steady state intensity shifts the mean cav- 
ity length. Introducing A^ = A^ — 2(7o/3, the intracav- 
ity amplitude is a^ ~ — i\/2K£i/(iA'j^ + k), where it has 




FIG. 1: A continuous variable pure state can be orthogonal- 
ized by coupling with a qubit via the Jaynes-Cummings (a) 
or the beam-splitter (b) interaction and then measurement 
of the qubit. Alternatively, simultaneously using the beam- 
splitter and two-mode-squeezing interactions can be used for 
state orthogonalization. This can be realized with cavity op- 
tomechanics to coherently manipulate the quantum state of 
motion of a mechanical oscillator (c). (PBS: polarizing beam 
splitter, FR; Faraday rotator). One of the drive fields is blue 
detuned and gives rise to a phonon-number-increasing process 
whereas the other is red detuned and gives rise to a phonon- 
number-reducing process. This is shown in (d), a truncated 
energy level diagram of the optomechanical system where the 
left kets describe the intracavity photon number and the right 
kets describe the mechanical phonon number. Each drive gen- 
erates a sideband at cavity resonance, which is shown in (e), 
an optomechanical spectrum. Thus, after erasure of the polar- 
ization information, photon detection at this frequency causes 
the mechanical element to undergo a coherent superposition 
of phonon addition and subtraction. 



been assumed that £ varies much slower than k. In the 
proceeding discussion this change to the detuning is ne- 
glected as the effect is small and can be readily com- 
pensated by frequency stabilization and/or appropriate 
pre-detuning. We turn now to the noise operators and for 
brevity solve the dynamics for a single drive frequency. 
We enter the mechanical and optical rotating frames via 
a — >■ ae~*^* and b -^ 6e~"^"', respectively. Assum- 
ing K <^ ujm, we make the rotating- wave approxima- 
tion and obtain a — iggyNab^'^' + \/2Kain — na and 
b = j(7oV^Q!^*'-'a'-'^\ where the brackets in the super- 
scripts are used to describe the two detunings we_consider 
(A = -I-WM7A = —(jJm) respectively. For go\/Na ^ k 



we use the adiabatic solution a ~ i^^/Nab^'^'^ + C, 
where C{t) = \/2^X!^dt' e-''(*-*')ai„(i'). The photon 
number scattered by the optomechanical interaction is 
n — /odiaQ^^aQ^^, which has been approximated to in- 
clude detection up to the drive duration t^k~^ and 
Oq^j — \/2Ka — fljii is the cavity output. For the h polar- 



ization driving the beam-splitter interaction (A = 
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fiV,|a,r and 
&Q is the mechanical field operator at the beginning 
of the interaction, time t — 0. For the v polarization 
(A = — wa/), which drives the two-mode-squeezing in- 
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We now consider 



weak drive such that the probability of more than one 
quanta being scattered is negligible. In this case, from 
the scattered photon number expectations, we intro- 
duce an effective beam-splitter parameter | = \f2G~hJ 
and an effective squeezing parameter r — ^JIG^t |28j and 
we describe the interaction using the effective unitary 



f/eff = 1 + (\ci\he~^^ 
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Here, 



and Lp arc the beam-splitter and two-mode-squeezer 
phases, respectively, which can be controlled via the 
phase of the drives. The fields at cavity resonance gen- 
erated via C/eff are spatially combined and filtered from 
the drive fields. Next, to control the weighting of iden- 
tity in the operation a weak displacement of amplitude 
[I is performed |3D|. Doing this to the h polarization. 



[/< 



cff 



1 + {\a\he-"^ - ra\h'^e"f -K 



^a\ 



this point the polarization of a scattered photon reveals 
how the phonon number changed. The field then passes 
through a wave-plate that performs ah — > ^(o/i + ^.v) 
and Oi, — >■ -jTiio-v — o.h) and is then incident upon a po- 
larizing beam splitter to conceal this information and 
allow for a quantum superposition. Conditioned on a 
h photon detection, the resulting mechanical state is 



T,Pi?/TT/Pr(/i), where Pr(/i) = TrMlT^T^p^) 



is the probability of photon detection and 
1 



T. 
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A V photon detection gives a measurement operator of the 
same form, however, with a n phase shift on the identity. 
Applications. — T^ provides a method to prepare and 
manipulate quantum coherence between the mechanical 
energy levels. Setting /x = 0, § = f and (j} = ip = -d + 7r/2 
we obtain the quantum state orthogonalizer T|_ — 
^(5g-.(^+'r/2) j^ 6te»(''+^/2))/y2 = rP|f . This quadra- 
ture is depicted in Fig. 2 as is its action on a displaced 
squeezed state. Such orthogonalization is heralded by the 
detection of a single photon that occurs with probabil- 
ity Pr(ft,) — 'I''^{{Pm Y), which is greater than zero for all 
physical states [3T] . We also note here that for states with 
zero phase-space mean, i.e. (i/'| h IV') = (V'l ^^ IV') — 0: one 
can also interpret these expressions as quanta subtrac- 
tion or addition to the state |V') yields a state which is 
orthogonal to 1-0) ■ Addition alone can orthogonalize all 
such states with a heralding probability of r'^{{b^b) -l-l)/2 
whereas subtraction alone has a heralding probability of 
{^)'^{b^b)/2. We thus further note that the operations 
6-/3 and b^ — P* can orthogonalize all pure states with 
(6) = /3. These operations may be simpler to experimen- 
tally implement, however, are less versatile as complete 
information of the state's mean is required as opposed to 
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FIG. 2: An equally weighted superposition of quanta ad- 
dition and subtraction can orthogonalize any pure quantum 
state, (a) The orthogonalizer T_i is a quadrature perpendic- 
ular to the angle i9 made by the input state's mean in phase 
space. The Wigner function (blue-cyan: positive, red-yellow: 
negative, larger ticks mark the origin and they increment by 
unity) of a displaced squeezed state (b), which has been or- 
thogonalized (c). A superposition of an initial state with an 
orthogonal state may be prepared to create a qubit from any 
initial pure state. In (d) such a superposition is shown by 
action with Tx -|- |m| e'^^'^^/V^, where |^| = r. 



the partial knowledge required by T,^. Returning to i2h 
one can now form a superposition of orthogonalization 
and identity, T^ — ii/^/2 + T±, to prepare a superposi- 
H.c). At tion of the initial state and an orthogonal state, i.e. a 
qubit, see Fig. 2(d). 

A mechanical resonator is a stationary system that al- 
lows T^ to be conveniently performed multiple times. 
Moreover, as the superposition weightings can be 
changed between applications this provides considerable 
versatility for quantum state engineering and quantum 
control protocols. For instance, one could realize the 
protocol by Dakna et al. [32] to synthesize an arbi- 
trary mechanical motional state. As another applica- 
tion, here we show that with N applications of T^, 
(2) one can transform the state \ij}) = ^^^ ■(/'„ \n) into any 

target state \(j)) — ^„ </>„ jn), i.e. arbitrary quantum 
state transformation. Our method uses only the sub- 
traction and identity components of T^ [53] and pro- 
ceeds in a manner similar to Ref. |32j and generalizes the 
scheme presented in Ref. [Mj. Specifically, by applying 

* = nf=i(Mj + ^3^)1^ = Eto C^b\ where v = f e"*-^, 
to the state \'ip) one can obtain \(j)) provided that the set of 
coefficients Ci is such that X]j=o" ^* i'i+n\/{i + n)\/n\ = 



Determining Gi can be readily achieved via ma- 
trix inversion and a solution exists provided that ^n ^ 
[29] . For the initial state having (6) = 0, the prob- 
ability of successful quantum state transformation is 

nti[(i)'(^^^)» + Im^|']/2, where {b%), is the phonon 
number expectation prior to the ith pulse. This prob- 
ability may seem low, however, the experiment can read- 
ily be performed with a megahertz repetition rate using 
a ^ 100 MHz mechanical oscillator and thus a practical 
number of heralding events can be attained in a reason- 
able time. Also, if the target state has a larger (smaller) 
dimension than the initial state one can apply creation 
(annihilation) as many times as necessary in order to 
make the dimensions the same prior to using $. 

An experimental approach. — There are numerous re- 
alizations of optomechanical systems and much progress 



has been made that can be built upon; the most per- 
tinent being Refs. [231 ISS] where phonon addition and 
subtraction were reahzed as separate operations. Com- 
bining these operations into a coherent superposition can 
be achieved with the setup in Fig. 1(c). Here we present 
an alternative route to fulfill the requirements of our pro- 
posed scheme using a mechanical element with a bulk 
acoustic wave vibration that forms an end mirror of a 
Fabry-Perot cavity (Mj. This configuration has the ad- 
vantage that the cavity decay can be controlled inde- 
pendently of the mechanical properties and such vibra- 
tional modes offer high mechanical resonance frequen- 
cies [37j. Moreover, simultaneous high reflectivity and 
high mechanical quality can be realized with multilayer 
crystalline reflectors 'SE]. A 40 /im diameter and sev- 
eral micrometer thick mirror has a mechanical resonance 
ujM /2t: = 200 MHz with a 20 ng effective mass. With 
a finesse of 5 x 10^, to achieve resolved sideband oper- 
ation, i.e. u!m/k= 10, a 75 /im cavity length can be 
used. For a drive laser with wavelength 1064 nm and 
a pulse duration of one hundred mechanical periods an 
optical power during the pulse of 1.3 mW is needed to 
achieve r^ =0.01. During the interaction the mechanical 
resonator also interacts with its thermal environment. To 
neglect the effects of environmental coupling we require 
that ^ = (fi/Q)(Ta;A//27r) <^ 1, where fi is the mechan- 
ical phonon occupation in thermal equilibrium and Q is 
the mechanical quality factor. For Q = 10^ and a 100 mK 
bath, which can be readily achieved using dilution refrig- 
eration, ^~ 10~^. Following the interaction the sideband 
needs to be separated from the drive field(s) prior to pho- 
ton detection. For higher mechanical frequencies the fil- 
tering requirements simplify. However, it is possible to 
achieve sufficient filtering even for a 200 MHz mechani- 
cal frequency using an optical displacement and spectral 
filtering (39j. Realizing the displacement with optical- 
fiber-based components, which provide excellent spatial 
mode matching, one can achieve an interferometric vis- 
ibility of 99.99% that suppresses the drive by 10**. The 
remaining drive can be further reduced by filtering with 
a cavity that has the same resonance frequency as the op- 
tomechanical cavity. To achieve a drive transmission 10^ 



times smaller than sideband transmission, a filter cavity 
amplitude decay rate of 2 kHz is required ^ . We would 
also like to emphasize that our scheme is robust against 
optical loss and inefficient detection as an optomechan- 
ically scattered photon that goes undetected does not 
trigger T^, hence, the primary effect of loss is to merely 
reduce the heralding probability [¥I] . To characterize the 
mechanical motional state, as the parameter regime con- 
sidered here is suited for the beam-splitter interaction, 
quantum state transfer of the mechanical motional state 
to the light 24, 42J can be performed followed by opti- 
cal homodyne tomography. This interaction, following 
action(s) with T^ to the stationary mechanical element, 
also provides a route to prepare optical continuous vari- 
able qubits or to synthesize arbitrary quantum states of 
a travelling optical field. 

Conclusions. — A superposition of quanta addition 
and subtraction can orthogonalize any pure continuous- 
variable quantum state with known angle made by the 
mean of the state's amplitude in phase-space. Such a 
superposition in combination with a controllable amount 
of the identity operation provides extensive control for 
quantum state engineering and quantum information ap- 
plications. For stationary systems it is convenient to ap- 
ply this tool multiple times, which we have utilized to 
illustrate how to perform arbitrary quantum state trans- 
formation. As the interactions we have used are available 
in many of the facets of quantum optics [33] , the tools we 
introduce can be realized in numerous physical systems. 
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Supplementary Material 



2. Arbitrary Quantum State Transformation 



1. Determining the Effective Unitary Interaction 

Central to our discussion in the main text is the mea- 
surement operator T^, which is used to describe the 
operation to the mechanical resonator via the optome- 
chanical interaction and then single photon detection. 
T^ ~ (1, 0| t/off |0, 0), where the ket is the initial state of 
light at the cavity resonance for the two orthogonal polar- 
izations used, the bra describes a /i-polarization photon 
detection with no v photon detection, and t/eff is the ef- 
fective optomechanical interaction including the manipu- 
lations to the optical field made after interaction with the 
mechanical resonator. In this supplementary we provide 
a discussion how [/off is obtained. 

The time evolutions described in Eq. (2) of the main 
text are generated by the beam-splitter and two-mode- 
squeezing effective interaction Hamiltonians. In the for- 
mer case a accumulates correlation with h and in the lat- 
ter case a accumulates correlation with 6^. For vacuum 
on the input of mode a, the expectation of the number 
operator in the output of mode a for the beam-splitter 
and two-mode-squeezing interactions are 
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(6^6), and sinhV(&fe^), 



respectively, where sin^(|) is the (intensity) reflectivity 
of the beam-splitter and r is the squeezing parameter. 
In the optomechanical scheme we have considered, the 
mean photon number scattered by the optomechanical 
interaction for the beam-splitter and two-mode-squeezing 
interactions are 
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){bX), andK) = (e2G.r_i)(5^fot), 



respectively. For small f , r, and Gt we then have 



f = ^2GhT, and r - y/2G^, 

for the effective optomechanical beam-splitter and two- 
mode-squeezing parameters, respectively. It is noted 
here that computing the mean number output in mode 
b can also be performed to yield these parameters. As 
both the beam-splitter and two-mode-squeezing pro- 
cesses are driven simultaneously, we expect that the 
effective optomechanical unitary take the form Ucs = 
exp[—j^{HBS + Hsq)t^, where Hbs oc a^b + ab^ and 
HsQ (X ab + a^b^ are the beam-splitter and two-mode- 
squeezing Hamiltonians respectively. To first order in 
the beam-splitter and squeezing parameters the effective 
unitary describing the cavity optomechanical interaction 
is then 

[/off = 1 + i^albe-''!' - ralb^e''^ - H.c). 

Finally, to obtain the effective unitary used for the mea- 
surement operator, the polarization manipulations to the 
optical fields, as discussed in the main text, must be per- 
formed. 



In the main text we introduced a scheme for arbitrary 
quantum state transformation that generates a target 
state from a known input state. Here we further dis- 
cuss our protocol and provide a specific quantum state 
transformation example. 

The protocol works as follows. For a known initial 
state 
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= E^" 



n=0 



which has no excitation beyond N quanta (or has been 
approximated by truncation at this level), any target 
state of the form 
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can be generated by applying a controllably weighted su- 
perposition of identity and subtraction N times, i.e. 
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<i' = l[{t^J+l^Jb)/V2 = J2c^b\ 

j=l i=0 

Applying this operation to the initial state we have 
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where we have used b \n) — ^/n \n — 1). 

The operation $ is a non-unitary process and the un- 
normalized matrix elements of the state after application 
of $ are 
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The target state |0) is reached when (n| $ |-0) — 0„. 
Provided that ^jv ?^ a set of coefficients Ci fulfilling 
{n\ $ 1-0) — (t>n can be determined via matrix inversion. 
Once a set of coefficients Ci is determined, a set of com- 
plex coefficients jjlj and Vj that satisfy (|3|) can also readily 
be determined via matrix inversion. 

We now provide a specific example of a quantum state 
transformation. Starting with an initial state \^) = |4) 
we wish to reach the target state \(j)) ~ (|1) + \4))/^/2. 
This target state can be reached with three applications 
of identity and subtraction. Solving Q we find that 
Co = \/24C3 and Ci = C2 = 0. As identity has been 
used with each application we set /i = 1 and obtain 

I'l + i^2 + t^3 = 0, 

1^1 1^2 + f^l^3 + ^2^3 = 0, 
VlV2Vj,V2l = 1. 



These equations can be readily solved exactly to pro- 
vide the relative amplitudes between identity and sub- 
traction to produce the target state. Numerical approxi- 



mations to the solutions and the intermediate states dur- 
ing the quantum state transformation process are shown 
in Fig. [31 




l^) = |4) 



= (|l) + |4»/v^ 



Ti = 1 - (0.294 + 0.510i)6 Ta = 1 - (0.294 - 0.510i)& 



T3 = 1 + 0.5896 



FIG. 3: An example quantum state transformation. Shown are Wigner functions (blue-cyan: positive, red-yellow: negative, 
larger ticks mark the origin and they increment by unity) of an initial Fock state (left) to a target state (right). The target 
state is reached by a sequence of three operations of a controUably weighted superposition of identity and subtraction. The 
relative amplitude between identity and subtraction for each step is shown. 
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